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Abstract 



We consider the models of relaxational dynamics within the framework of Berend- 
sen's and Nose — Hoover's thermostats. On studying the crucial case of ideal gas 
we come to the conclusion that both models mentioned above do not allow for 
describing the true thermodynamical equilibrium. 
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1. Introduction 



The thermodynamical concept of thermostat supposes that a system under con- 
sideration is a small part of a large system, or thermostat, which is big enough 
so as to be insensitive to changes inside the small one, and thus a very large 
number of particles is generally required to the effect. On the other hand, the 
small system may be in its turn rather tiny, even comprise only one particle. 
The thermostat's being macroscopical, creates serious difficulties for numerical 
modeling thermodynamical systems at fixed temperature and therefore, there is 
a need for simplified models that could accommodate the requirement of fixed 
temperature. The most straightforward approach to the problem is the use of 
the so-called collision thermostat, that provides adequate means for describing 
the system's interaction with ambient media. This method has its drawbacks 
in that it is time consuming and puts severe constraints on the precision level. 
Hence, it is very tentative to find certain simpler ways of modeling, using sys- 
tems that do not involve extremely large number of particles or stochastic equa- 
tions. The most familiar models in this respect are the Berendsen, Q and the 
Nose — Hoover, one. The idea of both is the same: we shall introduce in the 
right hand sides of a system under consideration relaxation terms that contain 
a dissipation parameter that should mimic the control exercised on the system 
by thermostat. To this end one uses a relation, of the most simple form pos- 
sible, that requires that the kinetic energy of the total system should be close 
to that prescribed by the heat energy determined by thermostat parameters. 
Thus, from the mathematical point of view, equations obtained in this way are 
an extension of the initial ones through the introduction of nonlinear dissipative 
terms in the right hand side. The prescription is quite general, appealing intu- 
itively, and easy to implement. But it reduces the problem to that of a finite 
nonlinear dissipative system, with no macroscopical features left and essential 
thermodynamical properties being muted. Cosequently, one may wonder as to 
whether the physical nature of the phenomenon has been preserved. 

In this paper we are going to consider examples of applying the prescription to a 
simple, but extremely important, case of ideal gas. Its analysis brings us to the 
conclusion that the Berendsen and Nose — Hoover systems allow for interesting 
nonlinear regimes, even for this simple case, but do not generate the system's 
thermodynamical equilibrium. 

2. Mechanical model 

Let us recall the construction of the thermostats mentioned above. We shall 
begin with the Berendsen one. Suppose we have a system of particles whose 
dynamics is governed by the Newton equations 
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= -—U(r) 
on 



The presence of a thermostat at temperature T is claimed to be accommodated 
through the introduction of non-linear dissipative terms in the equations indi- 
cated above, that is 



d 

niifi = -^—U(r) - jfi 
on 



(1) 



with the dissipative coefficient 7 being of the form 



7 = a 
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where N is the number of particles. Intuitively, it is "clear" that the model 
described above by equations ([!]), (^) should guarantee that the system be at 
temperature T, or at least its total kinetic energy be close to that value. The 
distribution of the kinetic energy per particle is not so clear. To be more precise 
equation (Q) in the notations of paper ^] contains the kinetic energy in the 
denominator; but it is not important for what follows. 

In the Nose — Hoover thermostat the dissipative variable 7 is treated on equal 
footing with other dynamical variables, and instead of the algebraic equation 
(^|) they write the differential equation 
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The system of equations (|T|) and (|^) is also quite appealing intuitively and looks 
excellently physically motivated. Again, at first sight it seems that the kinetic 
energy should be close to the heat one, at least if the system is kept long enough 
in the regime, and again one may have certain doubts, from the very beginning, 
as to its distribution with respect to degrees of freedom. 



2. Ideal gas 

Let us consider an ensemble of N nonintcracting identical particles of mass 

mi — m, i = 1, 2, . . . , N 



3 



that is an ideal gas, in a thermostat. Consider first the Berendsen one. Then 
the equation [l] takes the form 

mfi = —77^ (4) 

in which i takes the values 1, 2, . . . , N. It is more convenient to use the momenta 

Pi = mfi, i = 1,2,..., N 
and cast equations (|l|) in the form 

Pi = -JPi (5) 

On multiplying equations in (||) with pi,i — 1, 2, . . . , N and next taking their 
sum, we obtain the equation for the energy E 

dt \3k b TN J v 1 

In equilibrium the value of energy is fixed, and therefore, the right hand side 
of the equation given above be equal to zero, so that we obtain the two values: 
(1) E = 0, which is meaningless, and (2) E = ^k^TN , which is the right one. 
But the second requirement still remains, that is the existence of Maxwell's 
distribution for velocities of particles. Here we run across serious difficulties. 
The main point is that thermostat ought to bring the system to the equilibrium, 
even if its initial configuration is outside the latter. Let us consider an initial 
state of the gas that comprises two beams of particles, the first one comprising 
of -/Vi particles, all of them having the same momentum Pi , and the second one 
of the remaining N2 = N — N\ particles, all having the momentum P2, P\ 7^ 
P2. It is to be noted that equations (||) and (|^) are invariant with respect to 
permutations of the indices i, that is we may interchange the particles so that 
solutions to equations ([|), (|^), remain valid. This is a form of Gibbs' principle, 
that is particles of ideal gas being not distinguishable from each other. But the 
particles of the first beam have the same initial momenta, Pi, and therefore 
they are described by the same equation in which 

Pi=Px, i=l,2,...,JVi 

and the same is true for the particles of the second beam, in which 

Pi = P 2 , » = 2Vi + l,JVi + 2,...,JV 

Thus, we may write only two equations for the momenta Pi, P2 instead of N 
equations (H), that is 
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(7) 



From the last equation we infer that the locus of stationary states for a gas 
having an initial state consisting of two beams of N±, N 2 of particles, respectfully, 
is the ellipse in plane P\ , P 2 , given by the equation 



— Pf + f P| = 1 (8) 
3fc fc T7V 1 3mk b TN 2 v 7 

In Fig.l we illustrate the specific case of ten particles, the first beam comprises 
3 particle and the second one 7. The particles move on a straight line with the 
momenta P\ and P 2 , respectfully. The states given by equation^) are by no 
means equilibrium thermodynamical states. 

Let us turn to the ideal gas in the Nose — Hoover thermostat. On using 
momenta we may cast equations (Q), (|^) in the form 



Pi 

7 



(9) 
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Consider again the two beams of N\ and N 2 particles having initial momenta 
P\ and P 2 , respectfully. Similarly to the situation considered above, we may 
write equations (^|) in the form of three equations 
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Equations ([L0|) have stationary states given by the constraint 

7 = 

and equation(||). In this sense the Nose — Hoover thermostat is similar to 
the Berendsen one; they are also similar in that the stationary states are not 
thermodynamical equilibria, as is illustrated in Fig. 2. But it is worth noting 
that Nose — Hoover's thermostat, at least for the case of ideal gas we have 
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been studying, does not allow for attractor regimes, and in this sense it is even 
farther from the real thermostat than Berendsen's one. 



3. Conclusion 

The example of ideal gas shows quite convincingly that neither the Berendsen 
nor the Nose — Hoover thermostat can serve as a device for modeling ther- 
modynamical equilibrium. One may suggest, of course, that for system more 
sophisticated than an ideal gas they still could function to the effect, just owing 
to the complexity of a system under consideration, certain models of macro- 
molecules, for example. 

All the same, it is worthwhile to note that the Berendsen model isan interesting 
example of nonlinear dissipative problem, and its dynamics is characterized by 
the presence of a set in its phase space that attracts trajectories, that is an 
attractor. The latter is determined by the condition given by kinetic theory 



In fact, it has also a more subtle structure of preserving the beams of particles 
with equal momenta, as was discussed above. Thus, at equilibrium the states of 
the system lie on a sphere in phase space. If the equilibrium had been of thermo- 
dynamical nature, we should have had a distribution of momenta corresponding 
to the canonical ensemble at a temperature prescribed by the thermostat. In- 
stead, we have the equation (||), which describes the values of momenta for 
beams of particles of the same momentum. Thus, there is no canonical ensem- 
ble in the mechanical thermostat, which resembles more nonlinear dissipative 
systems like the Brusselator, M. At any rate, the system arrives at a stable 
state, following its evolution in time. 

In contrats, Nose — Hoover's system for ideal gas does not have attractors, even 
though the set of its stationatry solutions is stable (see Fig. 2 and 3). 
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FIGURE CAPTIONS 

Fig.l 

Ten one dimensional particles in the Bercndsen thermostat. The set of initial 
conditions for the momenta Pi , pj of the two beams of 3 and 7 particles corre- 
sponds to the set of trajectories which indicate the descent of the beams on the 
stationary states belonging to the ellipse given by equation (0). 

Fig.2 

Ten one dimensional particles in the Nose — Hoover thermostat. The horizontal 
axes are P%, P2, the vertical one gives values of 7. The closed curves indicate 
oscillatory motion round the set of stationary states, which is stable. The par- 
tition of the system into two beams is preserved, so that no thermodynamical 
equilibrium is achieved. 

Fig.3 

Ten one dimensional particles in the Nose — Hoover thermostat. The horizontal 
axes are Pi, P2, the vertical one gives values of 7. A large amplitude motion of 
the beam; no thermodynamical equilibrium. 
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Fig.3 




